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Abstract. We propose a scheme of ".R-type" photoassociative adiabatic passage 
(PAP) to create polar molecular condensates from two different species of ultracold 
atoms. Due to the presence of a quasi-coherent population trapping state in the scheme, 
it is possible to associate atoms into molecules with a low-power photoassociation 
(PA) laser. One remarkable advantage of our scheme is that a tunable atom-molecule 
coupling strength can be achieved by using a time-dependent PA field, which exhibits 
^) , larger flexibility than using a tunable magnetic field. In addition, our results show that 

the PA intensity required in the " -R-type" PAP could be greatly reduced compared to 
that in a conventional "A- type" one. 
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1. Introduction 

Recently, the realization of ultracold polar molecular gases has been regarded as one 
of the most promising research directions in the field of atomic and molecular physics 
[U E]. Ultracold polar molecules, with their long-range and anisotropic dipole-dipole 
interactions [3j H], have attracted much attention in a variety of research areas, such as 
quantum information science [5]- [7J and precision measurement [8]-|13j. 

There are two typical routes to achieve quantum degenerate gases of molecules. 
One is through the direct cooling of molecules, which is hard to achieve due to the 
complex internal levels of molecules [H]. The alternative one is to couple a pair 
of degenerate atoms by photoassociation (PA) [15] or Feshbach resonance (FR) [16J. 
However the diatomic molecule formed by a PA or FR process is usually loosely 
bound and energetically unstable. They have to be adiabatically transferred into a 
tightly bound ground state via a stimulated Raman adiabatic passage (STIRAP) [T7] . 
The success of the STIRAP is based on a coherent population trapping (CPT) state, 
which is accomplished by a pair of pulses in a counterintuitive sequence. During the 
adiabatic transfer, the system can follow the superposition between initial and final 
states, preventing any incoherent losses involving the middle unstable levels. Thereby, 
the high phase-space density of the initial gas can be coherently preserved. 

Currently, intensive experimental efforts to obtain quantum degenerate gases of 
molecules have been made by combining FR with STIRAP, which serves as an effective 
way to produce molecules in lower vibrational levels [18]- [21]. However, due to the 
strong vibrational quenching, a severe particle loss appears near the FR threshold. 
One way to solve this problem is to apply the optical lattice technique [HI H9) . in 
which inelastic collisions between molecules are well suppressed by preparing one single 
molecule per lattice site. Alternatively, all-optical transfer of molecules toward quantum 
degeneracy using a "two-color PA" method has been demonstrated experimentally [25] - 
[27] . where the excited molecules are moved down by a coherent dump field, instead 
of by spontaneous decay [28j-[30j. However, one common bottleneck with PAs is the 
small free-bound Franck-Condon factor (FCF), which requires an intense PA power to 
achieve an efficient adiabatic transfer [3TJ [32]. To date, the most promising way to 
overcome this PA weakness is by the FR-assisted PA scheme proposed first by Verhaar 
et. al. [331 El] and verified by many groups later on [35]- [39]. In terms of these 
studies, if a Feshbach quasi-bound state is adjusted close to the continuum, the atomic 
scattering wavefunction, acquiring some bound-state properties, becomes more localized. 
This gives rise to a dramatic enhancement of the free-bound FCF. As a result, the PA 
intensity required for a given atom-molecule transfer efficiency can be greatly reduced, 
compared with the case without the assistance of FR. 

In the present work, we propose an all-optical scheme to achieve a high transfer 
efficiency of atoms into molecules with a low PA power. For the purpose, we consider 
a " .R-type" atom-molecule conversion model (see figure [T] in solid arrows) through a 
photoassociative STIRAP procedure. Such a model is similar to the " i?-transfer" 
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suggested by Nikolov et. al. [40] as well as to the work by Band and Julienne [4"T] . 
In their works, molecules with an upper high-lying state are generated first through a 
step-wise PA excitation from free atoms (e.g. (Oi^) — > \m) — > \e) in figured]), followed 
by a radiative decay to populate a series of ground manifolds (e.g. |e) — > \g) in figure 
[T]). In this paper, we apply a coherent dump field to make the transition from state 
|e) to state \g), instead of by spontaneous emission. This leads to an accessible atom- 
molecule adiabatic passage between the initial (|0i j2 )) and target states (\g)). Compared 
with a conventional "A-type" model (see figure (TJ dash-dotted arrows), the CPT state 
supported in the "i?-type" scheme has been perturbed by a newly embedded state |m), 
which is absent in previous STIRAPs. As we will show, state \m) can help to reduce the 
power in PA field, whose stability properties will play a significant role in the molecular 
production. Under a simple numerical comparison, we have identified that the PA power 
required in the "A-type" model must be much higher than that in the ".R-type" model 
for achieving the same final efficiency 

This paper is organized as follows. In section 2, after briefly reviewing a similar 
idea of FR-induced STIRAP, we come up with our photoassociative STIRAP model and 
develop the underlying mean-field equations for the studies of a quasi-CPT description. 
In section 3, a generalized adiabatic theorem involving all the Bogoliubov collective 
modes is introduced to evaluate the adiabatic condition and quasi-CPT lifetime in 
our scheme. In section 4, numerical simulations for the cases described in section 2 
and 3 are implemented by using practical parameters. The laser profiles applied in 
the calculations are optimized according to adiabatic condition (section 4.1) and other 
relevant assumptions (section 4.2). Finally, a summary is given in section 5. 

2. Model and dark state theory 

Before moving to concrete illustrations of the model, we briefly review the idea of the 
FR-induced STIRAP method [121 S3]- In a typical FR-induced STIRAP, a number 
of colliding atoms undergo a strong association into quasi-bound molecules when a 
magnetic field is swept close to or across the FR, characterized by coupling strength a 
and binding energy e. Subsequently, these quasi-bound molecules are further transferred 
into stable molecules via a STIRAP. Clearly, here a and e correspond to QpA and S m 
(see figure [1]) in our scheme, respectively. A major advantage of our scheme is that 
flpA and S m can be manipulated more conveniently than a and e over the time scales. 
Because the latter quantities are highly dependent on the atomic intrinsic properties, 
especially the coupling strength a, which is fixed by the hyperfine interaction and is 
hence independent of time, e is experimentally tunable via an external magnetic field, 
while the former quantities, being controllable by optical means, are easily selected. 

Turning to our scheme, as depicted in figure [fl we study a "_R-type" five-level 
atom-polar-molecule formation. Two species of free atoms prepared in |0i) and | O2) 
states are first coupled into molecules in an intermediate high-lying state \m) with 
Rabi frequency flpA and detuning 5 m . Simultaneously, a pair of pump-dump lasers are 
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Figure 1. (color online) Schematic diagram: "iZ-type" transfer in solid arrows 
comprising the transitons of |0j 2) — ^ I™) ~ ^ l e ) ~ ^ |.<?) with the corresponding coupling 
fields CIpa, and fid, respectively; "A- type" transfer in dash-dotted arrows with 
|0i,2) — > |e) — > |<?) transitions characterized by Vt' PA and fi^- All the other parameters 
are described in the text. 



applied to move these loosely bound molecules in \m) down to the lowest molecular 
ground state \g), where fl p , Qd stand for coupling strengths and 5 e , 5 g for the two- 
and three-photon detunings. This scheme has several attractive properties. Firstly, the 
presence of state \m) brings one extra bound-bound transition from state \m) to state 
|e); hence it becomes easier for the PA field to associate atoms into the \m), rather 
than a higher |e) state. Secondly, transitions \m) — > |e) and |e) — > \g) are preferred 
because of the favorable bound-bound FCFs. Meanwhile, the (Oi^) — > \fn) transition is 
also accessible by an optimal control of the PA field in the time domain. 

As usual, we start our discussions with a set of coupled Gross-Pitaevskii's equations. 
In the mean-field treatment, where every quantum field operator has been replaced 
by its normalized amplitude ipi |44j . this yields 





2 V^m, 




itpo 2 = 


^PA ,* , 




tym = 




&PA , , 




- (S e + he) ~4>e ~ 


^ _ 9±d 
2 fm 2 ^91 


iip e = 


ilPg = 


-<w 9 - y^ 6 ' 





(1) 

(2) 

^e, (3) 

(4) 
(5) 

where ji (i=m,g) is introduced phenomenologically to describe the spontaneous decay 
of the \i) state to other undetected states, and it is possible to find a relatively stable 
\m) state with its decay rate 7 m <C 7 e [HI US]. The initial and target states are assumed 
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to be sufficiently stable with 7o 1(2) , g = 0. For an easy analysis without loss of the main 
physics, inter- and intra-species collisions have been ignored under typical parameters 
[16] . After a global gauge transformation, we can safely consider all the Rabi frequencies 
to be real positive values without loss of generality. 

A CPT state is always expected to move all the population into a target state as 
long as the adiabatic condition holds. In order to derive the corresponding adiabatic 
parameter, we first search for the CPT distributions for the following assumptions: 

V 0l , 2 = 0oe"^, ^ e = 0, W 9 = m , 9 e" 2 ^. (6) 

Here 0j is a steady-state amplitude, we consider </> 0l = 0o 2 = 0o for a balanced 
system, and /i is the atomic chemical potential. By ignoring all the decays and inserting 
equation ([6]) into equations ([I])-© with particle number conservation: 2(0 2 1 + 2 n -|-0 2 ,) = 
1, a generalized three-photon resonance is given by 

5 g ± = -2/i ± = , (7) 

Sm ± + Wa (3 + X 2 ) /2 

leading to the following CPT descriptions with e = 0: 



!> 2 

m 



:i+x 2 ), (8) 



= ( 9 ) 

i + ^i + nf M (3 + X 2)/2' 

09 = -X0m- (10) 

where flpA = tipA/S m , x — ti p /Qd- In equation (J7j), the choice of 5 g is determined 
by 8 m . If S m > 0, 5 g = 5 g+ and ji = fi + , whereas if 5 m < 0, 5 g = 5 g - and // = //_. 
From equations (18|)-( TT0|) . we note that when VLpa and \ both change from to large 
positive values, population initially prepared in states |0i j2 ) will be gradually converted 
into molecules in state \g) under three-photon resonance [equation (J7|)]. Also it is worth 
emphasizing that such a CPT state has been perturbed since m 7^ 0, and is called a 
"quasi-CPT" state. In the limit of VL P a 1, population in state \m) is virtually empty, 
we find that a complete transfer is still possible as long as x varies from to 00. In other 
words, the change by tip a has been accomplished by varying Q p ; thus the existence of 
state \m) is quite helpful for a relatively small tip a value. 

Actually, in the dynamics, if we use a strong PA laser to trigger the (Oi^) — > \m) 
transition, particle accumulations in state \m) will inevitably arise. Therefore, in order 
to avoid a considerable loss from state |m), pulse durations in STIRAP must be much 
shorter than \m) state's lifetime. On the other hand, if we deeply reduce the tip a value, 
the population in the \m) state will greatly be suppressed; meanwhile, a large fraction 
of atoms are left in the continuum, unpaired, because of a poor atom-molecule coupling 
strength. This conflict can be generalized to the properties of a quasi-CPT state, in 
which case one may prefer the use of moderate PA power. 
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Results in equations (|5])- (TTUI) are for the case of S m ^ 0. If 5 m = 0, i.e. the PA 
laser is exactly resonant with the free-bound transition, then equation ([7]) is reduced to 
5 g ± = ±tt PA /(6 + 2 X 2 ) 1/2 with the following CPT solutions: 

^ = 2^ = 2^/x 2 =(3 + x 2 ) _1 . (11) 

Equation fTTT]) shows a constant population ratio between states |0i i2 ) and \m), i.e. 
0o/0 m — 2. This equality contrasts with the standard CPT evolution, especially at t=0, 
which implies a poor transfer efficiency at 5 m = 0. As a result, a nonzero 5 m value is 
favored in our consideration. 



3. Adiabatic Theorem 



To derive the adiabatic parameter for the quasi-CPT state, we adopt a standard 
linearized approach as in jlTJ HH] by adding a small fluctuation 5ipi to the instantaneous 
steady-state solution 

V'Oi.a = (00 + Ho 1>2 ) Q (t) , V'e = tyeq 2 (t) , 1p m ,g = [<f> m ,g + &/w) q 2 (t) (12) 

where q (t) = exp[— / (i it') dt'], and (i(t) is a time- dependent chemical potential given 
by (i(t) = (see equation (JTj)). Substituting equation (f!2l) into the mean-field 

equations ([I])-© with the help of CPT descriptions, we eventually arrive at a set of 
linearized equations for the vector Sip = [Si/)q 17 Si[)q 2 , 6if) m , 6if) e , Si[)g\ T with its conjugate 
vector Sip* 
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the only nonzero elements. In addition, we have assumed detunings 5 e>g = — 2//(i) to be 
chirped [49]. 

Furthermore, we introduce a generalized Bogoliubov-de Gennes (BdG) equation for 
matrix M, 

M (t) Wi (t) = oJi (t) Wj (t) , (16) 

where Wj and Wj = [uj,Vj] T are the well-defined ith eigenenergy and eigenvector, 
respectively. Uj and Vj contain familiar Bogoliubov u-v parameters for each species 

u (v),. = [u (v) i0i , u (v) l02 , u {v) im , u {v) ie , u {v) ig ] T . (17) 

From the BdG equation, taking into account the special structure of matrix M, one 
can show the quantities u 2 are the eigenenergies of the matrix (A + B)(A — B), which 
can be obtained from the following equation: 

(utf ((c^ 2 ) 3 - a, (u; 2 ) 2 + a 2 (u 2 ) - a 3 ) = 0, (18) 

where the coefficients a« are given as 



Vt 2 + VL 2 

ai = (5 m + 2n) (5 m + 6/i) + \ \ (19) 



a-2 



2 



(2^i + Q 2 d (S m + 6/i)) + V P 16 dy , (20) 

a 3 = ( * m ffi " 3 (4^> + fig (5 m + 6/,)) . (21) 

Eigenenergies implied in equation (JTHjl comprise a doublet mode wo,i = and 
three pairs of excited modes (ujj, —Wjj (j=2,3,4). We find that ujj is real and has to be 
determined by biorthonormal relations for its corresponding eigenvector Wj. Detailed 
elucidations of biorthonormality have been published elsewhere [50] ■ In addition, we 
realize that the dynamical instability is impossible here due to the absence of collisions. 

To accomplish the goal of deriving the adiabatic theorem, we have to expand an 
arbitrary vector 5^/ in the dressed-state picture with a complete set of eigenvectors. 
By solving the BdG equation with eigenenergies, we are able to obtain w and Wi 
(unnormalized) explicitly using Gram-Schmidt orthogonalization, 



w = (-1,1, 0,0, 0,1, -1,0,0, Of, (22) 

( 00 00 - n -00 -00 , n \ T , OQ x 

V 2 0m 20 m 20 m 20 m J 

Here w , being a real dark state, is entirely decoupled with other eigenmodes, because 
its source term (wq 1 ^) and inter-coupling term (w^r^+w^) both vanish (see equation 
below for detailed notations), whereas, Wi is most likely to be triggered through 
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its nonzero inter-coupling strength (wj r^+w^), which can be roughly estimated by its 
decay rate 

< r »+« - ,wS:^^) s i/t * (24) 

where Q is a newly introduced vector complementary to wi with a well-defined 
normalization, 

wf?7 + Q = 1, (25) 

through the definition of 

MQ = wjv. (26) 

Here, v is a coefficient to be determined, and 77+ (and r/_ below) are given in 
[50J. Combining equation (1251 with ( 1261) . we find the vector Q takes a special form: 
Q = [%■> 9o> 9m j <?e, <?o ; 90; <?m; <Ze> <?g] T - Detailed expressions for q± and t> are presented 
in the appendix. 

The CPT lifetime r cpt defined in equation ( )24|) is clearly inversely-proportional to 
7 m and (f) m , which agrees with our intuitions. In other words, the presence of state \m) 
actually gives rise to a finite lifetime for the quasi-CPT state. Any pulse duration used 
in the system has to be much shorter than r cpt ; otherwise, a big particle loss from state 
I m) is unavoidable. One effective way to achieve a long T cpt is to search for a relatively 
stable \m) state with a small 7 m value. Other excited eigenenergies and eigenvectors are 
too complicated to list here, but they can be conveniently derived from equation (ITBl 
with (jig). 

Since other inter-coupling strengths for wi are also proportional to 7 m as in equation 
( 1241) and 7 m is considered to be much smaller than j e , we shall safely ignore the 
contributions from Wi and expand S^f in the parameter space with the help of three 
excited eigenmodes Wj (j =2,3,4) only, taking the form of 
4 

W = E (c^+w. " SVw*) . (27) 
i=2 

Through inserting equation (T27|) into ( JT3|) . and with the help of biorthonormality 
relations, finally, we obtain a set of coupling equations for Cj (£), 

+ iwj-c,- + wjr5* = -wj* (28) 

Terms like wj5^ have been eliminated in equation ( 128]) because the eigenvector 
Wj changes very slowly in the adiabatic limit. Generally speaking, if a system is said to 
operate in an adiabatic regime, population in any excited mode (nonzero eigen-mode) 
remains small. Hence, we shall introduce a typical adiabatic parameter definition 




r(t) = f^ ^ ^ «1 (29) 

A reduction in the r- value means an increase in the adiabaticity; in general, it can 
be accomplished by a longer pulse or a stronger pump field. In the adiabatic regime, if a 
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system evolves in a CPT state, an entire population conversion is achievable. However, 
the CPT lifetime implied in our model places a limitation for both the pulse duration 
and PA intensity, leading to a slightly larger r-value. This point will be discussed in 
section 4.1 toward the goal of obtaining an optimal pulse duration and PA intensity for 
an efficient transfer. 

In equation fl28|) . since dj can be ignored adiabatically, we further rewrite it as a 
series of linearized coupling equations: 

/ F + iuj G 
I G F - iuj 





(30) 
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up (j=2,3,4), 



fij (,'U'im'U'jm ^imVjrn) Te (Ui e Uj e Vi e Vj^) , (31) 

9ij 'Jm y^imMjm ^imV jra) Te (^ie^j'e ^ie^je) i (32) 

We solve c-, (c*) values from equations ( |30l) numerically and insert them into equation 
( |29i) . a time- dependent r-function is ultimately accessible. It needs to be noted that all 
the u(v)s in equation (}3~Tj) and (|32|) have been normalized according to biorthogonality, 

™Jv+™j = <%, w iV- w j = 0. (33) 



4. Numerical Analysis 

In the following numerical calculations, we intended to achieve a highly-efficient ground 
molecular production under an optimization of all the optical fields, including Qpa, 
Q p and Qj- From CPT descriptions [equations (IHlfTUI)]. we adopt a common pair of 
counterintuitive pump-dump pulses for \m)-\g) transition with the same width T 

^=^(l±tanh(^)) (34) 

where t Pjd are for the peak Rabi frequencies and central positions respectively. 

Based on equations f[5j)- f[TUj) . the PA Rabi frequency fipA, which must start from 0, is 
considered to share the same profile as Q p except for a different peak amplitude Q PA - 
Here, the detuning 5 m is fixed at a finite value for simplicity. 
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Figure 2. (a) r- value at t — t s versus pulse duration T under different PA amplitudes: 
from top to bottom Q° PA — 10 5 s _1 , 4 x 10 5 s" 1 , 10 6 s _1 and 10 7 s -1 , respectively, 
(b) r cp t estimated at t = t s as a function of ^p A - The four circles (from left to right) 
denote r cpt =17.8 ms, 2.2 ms, 0.97 ms, 257 /j,s with respect to the corresponding Slp A s 
shown in (a), (c) Final molecular production 77 = 2\ip g (oo)\ 2 versus pulse duration T 
under the same four cases as in (a). The other parameters are described in the text. 



4-1- Optimal PA pulse 

In what follows, we seek to gain from the r-value in equation fl29|) insights into the 
parameters, especially for an appropriate PA amplitude f2 PA and a pulse duration T. 
As we already understand, applying a longer pulse or a more intense PA laser will lead 
to a lower r-value. If a system's adiabaticity (r-value) is kept in a low level, which 
means the system will operate within the adiabatic regime, any excited modes are 
greatly suppressed. In a pure-CPT environment, adiabaticity indeed becomes a sufficient 
criterion for a complete transfer. However, in our scheme, we observe in addition to 
adiabaticity, a long CPT lifetime is another significant criterion for an efficient transfer. 

In a dynamical process, the r-value obtained from equation (1291) varies with time. 
We find that r-value estimated at t s which is defined by 4>l(t s ) = 2(f) 2 g {t s ) turns out to 
be a good estimate of the degree of adiabaticity. Thus, r s and r cpt values displayed in 
figure E(a), (b) are both evaluated at t = t s . 

Figure [2fa) and (c) present the variations of adiabaticity r(t s ) and final efficiency 
r](=\2ip g (oo)\ 2 ) as a function of pulse width T, respectively. As plotted in figure 2(a), 
either a longer pulse (from left to right) or a stronger PA amplitude (from the top to the 
bottom) leads to an improved adiabaticity. Furthermore, when fi Pj4 is very weak, such 
as 10 5 s _1 (in solid), the r-value is around 1.0, which cannot well satisfy the adiabatic 



11 



condition r < 1. Although at this time, the CPT lifetime in figure Wfo) is long enough 
(more than 15 ms) to support a longer pulse duration, a large part of atoms will be left 
in the continuum, resulting in poor molecular production, which is no more than 30% 
(see the solid curve in (c)). On the other hand, if we use an intense laser, fl° PA = 10 7 
s _1 , then the adiabaticity reduces into the 0.01 level, whereas simultaneously, the CPT 
lifetime is only around 250/is, leading to a dramatic reduction in t] as T increases (see the 
dash-dotted curve in figure 2(c)), because with a longer T value, a number of molecules 
decay spontaneously due to j m . Obviously, if T < 20/is, a relatively higher r\ value 
(>50%) is still attainable. 

In addition, we study two moderate cases with the PA amplitudes: flp A = 4 x 10 5 
s _1 (in dashed) and flp A = 10 6 s _1 (in dotted). No impressive differences are observable 
in adiabaticity according to figure [2(a), where both are around the 0.1 level. Meanwhile, 
the T cpt values represented in figure [2](b) are both close to 1 ms, which do offer more 
space for a tunable T value. Final efficiencies in figure [2]^c) clearly exhibit a T-dependent 
feature, while staying at a highly efficient level compared with two former cases. 

In light of the above discussions, we conclude that the adiabaticity indeed serves as 
a useful tool to select favorable parameters. Meanwhile, it is equivalently important 
to take the CPT lifetime into consideration. Here, we prefer to use T = 30 fis, 
n° PA = Ax 10 5 s _1 . 

4-2. Optimal pump-dump pulse sequence 

The goal of this subsection is to design the optimal pump-dump two-pulse sequence to 
maximize the yield of molecules. Clearly, there are five individual parameters to be 
determined: t p , t^, 0°, fl d and 8 m . Such a five-parameter variation is difficult to carry 
out. However, from the CPT descriptions, we guess that the population dynamics are 
most likely to be affected by the ratio % = ^p/^rf instead of the fl p and fid values. 
Therefore, we introduce two new variables, which are dt = — t p for pulse delay and 
X° = Qp/Qd, f° r peak amplitude ratio. In combination with the one-photon detuning 
5 m , there are three effective quantities to be optimized. 

Figure [3]Ja) and (b) show the contour plots of the final molecular productions for 
sets of [5 m ,x°] and [6 m ,dt], respectively, where the lighter areas correspond to higher 
efficiencies. Especially, pure white regimes are for rj > 80%. These two mappings 
have several attractive features. Firstly, a symmetric pattern along the 5 m direction is 
explicitly observable, which can be ascribed to the existence of a three-photon resonance 
[equation (JJJ)]. Regardless of whether 5 m is positive or negative, either 5 g+ = — 2/x+ or 
5 g - = — 2/i_ will be satisfied. In other words, a double-resonant condition must hold 
on both sides of 8 m , leading to a symmetric double-peak pattern. Similar patterns have 
been demonstrated by the Autler-Townes splitting effect [511 E2] , which usually takes 
place when an optical field is detuned close to an exact transition frequency. To be more 
understandable, if we artificially add a small perturbation to a resonance, the double- 
peak profile will be correspondingly shifted. Since this shift employs no improvement 
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Figure 3. Contour plot of final transfer efficiencies under different sets of (a) \S m , x°] 
where dt — 105/is and (b) [S m ,dt], where x° — 6. Here, 0° = fO 7 s -1 , t p = 75 fis. 
Lighter areas correspond to high efficiencies. Pure white areas denote efficiencies more 
than 80%. The other parameters are listed in the text. 



in the molecule production, we will leave this point for future interested readers. 

Secondly, if we fix \S m \ around 0.05 x 10 7 s _1 and gradually increase the values of x° 
and dt, the final transfer efficiencies express similar variations. Seen in figure |3]^a), based 
on = 10 s -1 , dt = 105/^s, if x° changes from 1.0 to 6.0, a dramatic enhancement for 
1] is explicit. When further increasing x° U P to 10.0, rj values will be slowly decreasing. 
A similar trend with rj as the pulse delay dt varies is depicted in figure [3](b) where 
X° = 6, = 10 7 s _1 . When 5 m = 0, efficiencies are very poor, which agrees with our 
CPT predictions equation (TTTj) in section 2. Finally, we find that the base value of 
offers few contributions to the transfer. If is set as 2 x 10 7 s~ 1 , we will obtain a much 
analogous contour plot to figure |3]^ a) (not shown). 

A brief conclusion for the sections 4.1 and 4.2 is that we are provided with rich 
ways to select relevant parameters for optimal atom-molecule conversion. 

4-3. Population dynamics 

In the following, we consider a concrete example in our "i?-type" scheme using the 
parameters based on our previous discussions. Optimal parameters are given by T = 30 
/is, tt° PA = 4 x 10 5 s-\ llj = 6x 10 7 s _1 , tt° d = 10 7 s -1 , t p = 75 /is, t d = 180 fis, 
\S m \ = 5.4 x 10 5 s -1 , 7 m = 3 x 10 4 s _1 and 7 e = 10 7 s _1 . Numerical results are plotted in 
figure HI By directly integrating the mean-field dynamic equations (JTJ)- (JSJ) , we produce a 
population dynamics which contains all the field amplitudes in figure H^b). Observably, 
more than 85% of the atoms ipo 12 convert into ground-state molecules ipg. Compared 
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Figure 4. (color online) (a) CPT dynamics, (b) population dynamics, (c) time- 
dependent adiabaticity values, (d) three excited eigen-energies as time changes where 
W3 is at least two orders smaller than W2,4 (shown in the inset). 



with the CPT dynamics shown in figure Ufa), a good agreement is clearly seen, except for 
a slightly lower 2ipg coming from spontaneous decays. In particular, we need to mention 
that the ^J^m) amplitude (dotted) has been deeply suppressed, with a maximum value 
smaller than 0.02. 

Figure H](c) represents the adiabaticity defined in equation (129!) as time changes 
(in solid). By solving equations (!30|) numerically, a complete r-value is able to be 
determined from the Cj values. Three excited eigenenergies obtained from equation (1181) 
are displayed in figure H](d) and the inset, where u 3 is smaller than o; 2i 4 by orders of 
magnitude. In the dressed state picture, Uj stands for the energy of the jth eigenstate, 
and generally speaking, a higher-energy eigenstate is usually more difficult to populate 
than a lower one. Thereby, in deriving the adiabaticity, we shall safely neglect the 
contributions from U2,4 and w 2i 4, simplifying equations ( )30l) with uj 3 , w 3 only, yielding 

iuj 3 c 3 + / 33 c 3 = -W3 (35) 

which leads to a reduced assessment on adiabaticity: r a = wjQ>/(iu) 3 + / 33 ) / y/3. 
Clearly, r a matches with r in figure 0](c) perfectly. 

One critical concern in our scheme is the stability of state \m), which indeed plays a 
vital role in determining the final transfer efficiency. In our calculations, we use the \m) 
state lifetime to be r m = l/{2n r ) m ) = 5.3/is, which is comparable with the earlier work 
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Figure 5. The final transfer efficiency r\ as a function of the intermediate state lifetime 
r m . An arrow corresponds to the lifetime we have used. 

of Napolitano et al [45]. Figure [5] displays how the final efficiency rj(=2if>g(oo)) varies 
as a function of r m . Clearly, if r m is smaller than lfis, rj drops rapidly as r m becomes 
shorter. However, if we are able to find a more stable intermediate state, with a lifetime 
longer than 10 fis, the corresponding transfer efficiency reaches as high as 90%. The 
arrow shown in figure [5] points to the r m value used in our paper. 

Finally, it is meaningful for a numeric estimate of the feasibility of our scheme 
by taking the KRb molecule as a possible candidate in experiment. Based on the 
predictions in [SUES], it is experimentally possible to find out an appropriate \m) state 
with a relative long lifetime, e.g. 3 1 S + . There is a transition dipole moment of ~ 1 eao 
for the state 3 1 S + , which corresponds to the lifetime of several /is. In addition, the l l U 
is a good candidate for the high-lying |e) state because of its purely singlet character and 
favorable transition dipole moment associated with the lowest singlet state 1 1 E + (\g) 
state). A rough estimation of the PA power, adopting the parameters of the free-bound 
FCF Ifcf ~10 -14 m 3//2 for KRb [5j]-[56] and an initial atomic density n = 10 20 m~ 3 , 
gives rise to a PA laser intensity of I pa = 2ceo{hQ°p A V> / fi) 2 ~ 512.7 W/cm 2 for our 
" R-type" scheme, where c is the light velocity, e$ is vacuum permittivity, \x is the dipole 
moment, Vp PA is the electronic Rabi frequency defined by £lp A = Qp\/ {lFCF\fno) [54"] . 
However, in a standard "A-type" system (see figured]), the absence of state \m) requires 
a more intense PA field to stimulate particles into the highly excited state \e). As a 
result, to achieve the same production rate of molecules as in our "R-type" case, the 
required PA power must be f2p A ~ 6.4 x 10 6 s -1 . The corresponding PA intensity is 
I' PA ~ 1.31 x 10 5 W/cm 2 , giving other parameters the same as in the "R-type" case. 
Evidently, the above numeric estimate shows a more than 250 times power reduction in 
our "R-type" approach. 
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5. Summary 

Although the magnetic FR-assisted PA technique has been considered as the most 
promising way for the purpose of overcoming the PA weakness, the primary drawback 
in such a scheme is the strong inelastic-collisional loss of Feshbach molecules, especially 
when the magnetic field is tuned near the resonant point. To eliminate the bottleneck in 
the magnetic FR-assisted PA technique, we work out a robust all-optical atom-molecule 
conversion model through a "i?-type" photoassociative STIRAP, where an intermediate 
state \m) is introduced to form a quasi-CPT state. In terms of the detailed adiabatic 
theorem, we show that the quasi-CPT state can lead to a higher atom-molecule transfer 
efficiency with a lower PA laser power, compared to the normal CPT state in the 
conventional all-optical "A- type" two-color PA configuration. 

The key reason for the lowered power of PA laser is due to the existence of an 
intermediate state \m). In this case, it is easier to photoassociate free atoms into this 
low-lying \m) state, rather than a high- lying |e) state. Since molecules in state \m) are 
unstable, the subsequent STIRAP transfer from the intermediate state \m) to the final 
state \g) must be rapid enough to avoid the loss of molecules from the \m) state, here 
characterized by a finite CPT lifetime. In addition, we also show that the lifetime of 
state \m) will significantly affect the final transfer efficiency A specific estimation to 
illustrate the feasibility of our approach is performed. Finally we want to emphasize that 
the scheme proposed here is the first one to overcome the inefficiency of PA with only all- 
optical fields involved. This may open up new opportunities for experimental endeavors 
to create polar molecular condensates directly from ultracold atoms. A more careful 
treatment taking into account nonlinear collisions will be left for future explorations. 



Acknowledgments 

This work is supported by the National Natural Science Foundation of China under 
Grant No. 10588402, the National Basic Research Program of China (973 Program) 
under Grant No. 2006CB921104, the Program of Shanghai Subject Chief Scientist 
under Grant No. 08XD14017, and the Shanghai Leading Academic Discipline Project 
under Grant No. B480 (W.Z.), the National Natural Science Foundation of China under 
Grant No. 10974057 and No. 10874045, Shanghai Pujiang Program under Grant No. 
08PJ1405000 (L.Z.). 



Appendix A. Q vector 

qi values in vector Q are given by 



go = 
q m = 



<j>l + 2<j>l (1 + x 2 ) 



2vVt PA <f)Q<f): 
1 



PAS 



(A.l) 
(A.2) 
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and 



* = H (A ' 3) 

% = -r^— (A.4) 
V\l PA <p m 



?§+ m± +x 2 ) 



Qpa4>: 



rn 
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